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ABSTRACT 

Together with their fingerprint modes, molecules carry coherent vibrations of all their 
atoms (phonons). Phonon spectra extend from ~20 to more than 10^/im, depending 
on molecular size. These spectra are discrete but large assemblies of molecules of 
the same family, differing only by minor structural details, will produce continua. As 
such assemblies are expected to exist in regions where dust accumulates, they are 
bound to contribute to the observed continua underlying the Unidentified Infrared 
Bands and the 21-/im band of planetary nebulae as well as to the diffuse galactic 
emission surveyed by the Planck astronomical satellite and other means. The purpose 
of this work is to determine, for carbon-rich molecules, the intensity of such continua 
and their extent into the millimetric range, and to evaluate their detectability in 
this range. The rules governing the spectral distributions of phonons are derived and 
shown to differ from those which obtain in the solid state. Their application allow the 
extinction cross-section per H atom, and its maximum wavelength, to be determined 
as a function of molecular size and dimensionality. Chemical modeling of more than 
15 large molecules illustrate these results. It is found that the maximum phonon 
wavelength of a 2D structure increases roughly as the square of its larger dimension. 
Spectral energy distributions were computed as far as 4000 /im, for molecules up to 
50 A in length. 
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1 INTRODUCTION 

Phonons in extended media have long been the subject of 
detailed investigations, in particular of their DOM (density 
of modes, g) and velocity, U, as a function of angular fre¬ 
quency uj (e.g. see Kittel 2005). Graphite, in particular, was 
abundantly investigated from the early work of Nicklow, 
Wakabayashi and Smith (|1972l) to the recent one by Mohr 
et ah (|2007p . By contrast, the study of vibrations in finite 
structures like molecules, has largely been limited to their 
fingerprints, i.e. those vibrations which involve only a hand¬ 
ful of closely bonded atoms within the structure (functional 
groups) and usually fall above 500 cm“^ (e.g. see Colthup et 
al. 1990). Molecular phonons differ from fingerprints in that 
they involve the whole structure; their spectra are quasi- 
periodic and extend from ~20 to more than 10^ jam (de¬ 
pending on molecular size) and so may contribute to the 
continuum that is observed to underly the UIB spectra of 
galaxies and planetary nebulae: see, for instance, Smith et 
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al. (|2007p . Compiegne et al. (|201ip . Zhang and Kwok (|201ip . 
Fischer et al. (|2014p . 

Low-frequency phonons are presently of particular in¬ 
terest to astrophysics as electromagnetic emission in the 
microwave/millimeter ranges has now been detected, mea¬ 
sured and surveyed all over the sky (e.g. see Finkbeiner et 
al. (|1999p : Ghosh et al. f (|2Q12p and bibliography therein). 
More recently, Ade et al. (|2Q14I) and Abergel et al. (|2014l) . 
analyzing data from the Planck astronomical satellite, ob¬ 
tained quasi-continuous spectra from 20 to 3000 GHz (100 
to 15000 jam). In this range, several emission sources appear 
to be in action: thermal dust, spinning and magnetic dust 
(Draine and Lazarian (1998), Lazarian and Draine (2000)), 
free-free electrons, synchrotron. Much effort has gone into 
combining these various sources so as to fit observations. The 
purpose, here, is more limited and focused on molecular vi¬ 
brations that may be excited either thermally or otherwise. 
We seek to define the physical (structure, size) and chemi¬ 
cal (composition) conditions to be met by a molecule for its 
low-frequency vibrations to contribute significantly to the 
observed spectra, and to determine the nature of the vibra¬ 
tions involved. This also differs from Papoular (|2Q14ap and 
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()2014b|) . in which the maximum wavelength was 1000 fim 
and only a mixture of molecules was considered but not their 
respective contributions to the emission. Also, this work can 
be considered a first step in extrapolating the extinction 
curve of carbonaceous dust into the microwave range. 

Fortunately, the study of molecular phonons is now 
made easier by computer modeling of finite, but very large, 
bonded structures (hundreds of atoms). In the course of such 
an investigation, it turned out that the frequency depen¬ 
dences of g and V, for the transverse acoustic phonons of 
ID and 2D molecules, differed from the well-known Debye 
laws for the case of bulk media (see Kittel 2005: p oc 
and V constant, in the limit of low frequencies, in which 
we are primarily interested here). While 2D molecules could 
be likened superficially to very thin plates, the theoretical 
treatment of the latter involves Young’s modulus and Pois¬ 
son’s ratio (see Reddy 2007) so it cannot be extended to 
molecules. 

Based on the study of the transverse and longitudinal 
acoustic phonons (TA and LA, respectively) of simple linear 
structures, it is shown here that the reason for this behavior 
of untethered molecules is that the vibrational mechanics is 
dominated by the force which resists bond bending, so the 
sound speed increases with frequency, by contrast with the 
constant sound speed in crystals (Sec. 2). 

In most cases, the topological complexity of the struc¬ 
ture does not allow direct analysis of motion for each in¬ 
dividual phonons; only the DOM is readily available from 
modeling. But what is needed to predict the lowest natural 
frequency of a given structure as a function of its size, will 
turn out to be the sound velocity. A mathematical analysis 
of the general case led us to a simple formula by means of 
which V can be derived from the sole knowledge of g. At 
low frequencies, g and V are governed by power laws; an¬ 
other outcome of our analysis is a set of relations between 
the dimensionality of the structure, the type of forces in ac¬ 
tion and the indexes of the two power laws which govern g 
and y, respectively (Sec. 3). Note that the mentioned di¬ 
mensionality refers here to the carbonaceous backbone of a 
molecule. 

In Sec. 4, the theory developed in Sec. 3 is applied to 
several 2D molecules which are shown to obey the predicted 
relations even though they may be bent and/or warped. The 
lowest phonon frequency (always that of a TA phonon) is 
shown to scale like the square of the major dimension of 
the molecule. In Sec. 5, the theory is further illustrated with 
several individual molecules. 

In Sec. 6, the integrated band intensity (absorption) 
spectrum is computed for several large molecules. The cor¬ 
responding extinction cross-section per H atom in the ISM 
(Interstellar Medium) is computed. The absorption spec¬ 
trum is used in Sec. 7 to obtain the emission spectrum of 
a particle in thermal equilibrium with its environment. The 
case where the particle is too small to reach thermal equi¬ 
librium is then considered, detailing a procedure used to 
obtain the corresponding emission spectrum, which is dis¬ 
tinctly different than in the previous case. It is concluded 
that reasonable amounts of C atoms locked in similar struc¬ 
tures are sufficient to emit detectable radiation as far as 4000 
A at least. 


2 PHONONS IN LINEAR ATOMIC CHAINS 

In the course of this work, a large number of carbona¬ 
ceous molecules were simulated by means of various algo¬ 
rithms of computational organic chemistry as embodied in 
the Hyperchem software provided by Hypercube, Inc., and 
described in detail in their publication HC50-00-03-00. The 
semi-empirical (mixed quantum/classical) PM3 method, ini¬ 
tially developed by J. P. Stewart ()I990D . was used through¬ 
out, as its algorithms are specially optimized for carbona¬ 
ceous structures; its accuracy for vibrational frequencies is 
a few percent. It does not require the initial selection and 
setting of particular force fields. 

The procedure for each molecule is as follows. First, 
the molecule is schematically represented on the computer 
screen. Then, the code is prompted to optimize the molecu¬ 
lar configuration by adjusting the various bond lengths and 
angles so as to minimize the total bond energy. At that point, 
a Normal Mode analysis can be automatically performed, 
which delivers the wave number and integrated band in¬ 
tensity (oscillator strength) of each vibrational mode of the 
system. 

As an example, consider a linear chain of 20 doubly- 
bonded carbon atoms, =C=. In order to impose a strict 
dimensionality of I, it is terminated at each end with =C- 
H; this is an ansatz provided by the simulation code under 
“Allow arbitrary valence” or “Allow ions”, and boils down to 
adding a local electric charge. Two alternative terminations 
were also used: CH2 and C=0, and gave similar results. 

The atomic motions were monitored on the screen for 
each normal mode found by the code. Three types of mo¬ 
tion were found: a) transverse atomic displacements parallel 
to the screen plane; b) transverse displacements orthogonal 
to that plane; c) longitudinal displacements parallel to the 
central line. The first two types have essentially identical 
frequencies, the third has distinct frequencies shifted higher 
up. In all 3 cases, the amplitudes of atomic excursions dis¬ 
play, all along the chain, a clear stationary sinusoidal pat¬ 
tern of standing waves which can be associated with har¬ 
monic waves traveling in opposite directions, whose wave¬ 
length along the chain decreases with increasing frequency. 
At each transition from one normal mode to the next, the 
number of half sinusoids abruptly increases by one. This 
is in line with the classical analysis of regular crystal lat¬ 
tices (see Kittel 2005). It is therefore possible to introduce 
here, too, the notions of phonon and wave vector k = I'k 
where \ is the wavelength of the phonon at the correspond¬ 
ing angular frequency, uj — 27r z/ c where v is the frequency 
(wavenumber) and c is the velocity of light in vacuum; the 
density of modes, ^(cc), like the density of states (DOS), can 
also be defined as the number of normal modes per unit fre¬ 
quency. By definition, the phase velocity is then V — duj/dk 
. If the normal modes are labeled by integer i, starting from 
I, in order of increasing frequency, then Ap(z) = 2Lmoi/L 
where Lmoi is the length of the chain; hence k{i) — ZTr/Lmoi- 

A careful monitoring of the vibrational patterns shows 
that the expression for Ap(z) must be divided by a factor /, 
which depends slightly on the molecular length and termina¬ 
tions, but remains in the range of 1.5 to 2 (here, 1.73). This 
is due to the fact that, in the present cases, the boundary 
conditions are not forced upon the molecule by stretching it 
between two fixed points, but are determined by the electri- 
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Figure 1. Transverse phonons. From bottom up: A) wave vector 
dots: from chemical modeling; line: fit with eq. 4. B) DOM 
g{y): crosses and pluses for parallel and perpendicular vibrations, 
respectively. C) V{u)\ crosses and pluses for parallel and perpen¬ 
dicular vibrations, respectively; circles: fit with eq. 12. All from 
the modeling of the chain of 20 doubly-bonded C atoms. 



Figure 2. Same as graph C in Fig.[T]for 7 linear chains of different 
lengths, including that of Fig. [T] (circles), except here, for clarity, 
only one of the two sets of the transverse vibrations is represented. 
Note the logarithmic scales, the wave numbers in abscissa and the 
gap near 440 cm“^. The rising part of the graph follows roughly 
the law V = 3.210^ . 


cal properties of the H termination at each end of the chain; 
as a result, the chain ends do not coincide with nodes of the 
standing wave. 

2.1 Transverse vibrations 

Figure [T]shows, as a function of the frequency, v: A) the wave 
vector, k] B)the DOM, g; C) the sound velocity, V = dcoldk, 
of the transverse vibrations, for atomic displacements paral¬ 
lel (x) and perpendicular (+) to the screen plane, all from 
the modeling of the chain of 20 doubly-bonded C atoms. As 
expected, there is a near coincidence of the 2 sets of points 
in each of graphs B and C . These vibrations, of course, 
are reminiscent of plate vibrations, and correspond to the 
transverse acoustic (TA) waves of solid state physics. 

Note the sound speed variation at low frequencies. 
This trend is confirmed by the modeling of six other lin¬ 
ear doubly-bonded carbon chains, made of 12, 20, 40, 80, 
120 and 160 C atoms, and terminated with CH2 at each 
end. The same procedure as above yields the TA sound ve¬ 
locities drawn in Fig. 2, together with that of Fig. [1] (for the 
sake of clarity, only one of the two sets of transverse vibra¬ 
tions is included). To the accuracy of the procedure, all 7 
sound velocities are nearly coincident, except for the chain 
of Fig. [H which is only slightly different, this being due to 
the different terminations. As expected, V is nearly indepen¬ 
dent of the chain length. The logarithmic scales highlight the 
power-law dependence of the sound speed on frequency at 
low frequencies, V = 3.2 10^ . 

Returning to Fig. [U note that the DOM goes to inhn- 
ity at = 0 and ~ 700 cm“^; cor relatively, by the dehni- 
tions of V and g, V goes to 0. These are the only singulari¬ 
ties predicted by Van Hove (|1953p for inhnite linear chains. 
With hnite chains, however, another type of singularity oc¬ 
curs near the minimum of the DOM at oj/2, characterized 
by a gap in mode frequencies, wider or narrower depending 
on the type of terminations. This is due to a chance reso¬ 
nance between one TA mode and any other strong mode of 


the system, which gives rise to a strong mixed mode. The 
latter then perturbs the distribution of neighboring modes 
by pushing them aside. This in turn, locally disturbs the 
DOM and sound speed, ah the more where the DOM is low, 
as is apparent in Fig. [1] In the present case, this occurs for 
the second longitudinal mode at 439 cm“^, which resonates 
with transverse modes at 442.6 and 443 cm“^, as shown by 
the following study of longitudinal acoustic (LA) branch of 
vibrations. 


2.2 Longitudinal vibrations 

In monitoring, on the computer screen, atomic displace¬ 
ments relative to the chain backbone, in each normal 
mode, it is straightforward to distinguish longitudinal from 
transversal vibrations. Fig. [3] shows the behavior of wave 
vector, DOM and sound velocity for the longitudinal vibra¬ 
tions of the same structure: by contrast with the case of 
transverse vibrations, the latter two are nearly constant at 
low frequencies. Note that, in the familiar case of 3D media, 
while V behaves likewise, g(p'), for its part, scales like 
instead. Remarkably, the sound speed of longitudinal waves 
is an order of magnitude higher than that of the transversal 
ones. As a consequence, the hrst longitudinal mode occurs 
at a much higher frequency than its transverse counterpart. 
The second is the one that causes the frequency gap in Fig. 
[U as explained above. There is also a gap in the longitudinal 
spectrum, but again at a much higher frequency, 855 cm“^. 
In this case, it is due to interaction with a strong CH bend¬ 
ing mode, which is clearly revealed on the modeling screen, 
at the chain end. 

In this work, modes involving strong motions of termi¬ 
nal heteroatoms, and extending only slightly into the chain, 
were observed on ah structures. They occur at the high fre¬ 
quency end of the transverse phonon spectral range, where 
end modes are strongly coupled to the last transverse mode. 
They are reminiscent of the “surface” modes studied by 
Wallis (|1957p . who showed that they can affect the thermal 
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Figure 3. Longitudinal phonons. From bottom up: A)wave vector 
dots: from chemical modeling; line: fit with eq. 6. B)DOM 
C)sound speed dots: chemical model; +: fit with eq. 

12. 


properties of very small powder grains. Here, however, these 
modes are relatively very few, so they were not given partic¬ 
ular attention, especially considering that we are primarily 
interested in larger grains and lower frequencies. 

Also note the absence of an optical phonon branch, as 
the atoms on the chain are identical. 


2.3 Vibrational mechanics of phonons 

Now, why should the velocity of transverse waves depend 
on the wave number, as documented above? A plausible 
surmise is that the main stress acting on the transverse vi¬ 
brations considered above is not the internal longitudinal 
tension (since the string is not stretched) but that which 
is induced by the bending of adjacent atomic bonds. It is 
known, indeed, that, for each atom type, there is a pre¬ 
ferred angular distribution of the valence electrons. In the 
present case, the double bonds on both sides of an atom re¬ 
sist any stress that tends to bend their alignment. Molecular 
mechanics models (CHARMM, Gaussian, etc.) account for 
this by introducing explicitly potential term in 0^, where 0 
is the bending angle. In semi-empirical modeling codes, like 
PM3 here, the effect automatically results from the quantal 
description of the valence electrons. In the present instance, 
by noting the increase, A t/, in potential energy upon bend¬ 
ing one of the pairs of bonds step by step from 0 to 5 deg, 
it was found that AU = eO^, where e = 1.85 10“^^ erg/rad^ 
(810“^ kcal/mol/deg^) and 0 is in radians. 

Let us assume that this is the only potential at work, 
and postulate for 0 a trial solution in the form of a traveling 
wave, expi(ct;t + rka), where r is an integer designating the 
rank of an atom in the chain, and a the distance between 
adjacent atoms, 1.28 A . Also, let yr be the transverse dis¬ 
placement of the rth carbon atom along the chain. Then, the 
angular displacement of the adjacent C=C bond relative to 
the chain’s backbone is 

ar — sina = - - - = sm [ka/2) -, (1) 

and 0 — 2a. The total potential and kinetic energies of the 
system are, respectively. 


1/ = 4e a^ and T = 


myr 

2 


( 2 ) 


where m is the atomic mass, and it was assumed that U is 
not dependent on the Ar’s. One can then write an indepen¬ 
dent Lagrange’s equation of the form 


d ST ST ^ SU 
dt S Otr S ar S ar 


(3) 


for each C atom in the chain, which gives the condition of 
existence of the trial solution: 


= 32 esin^(^). (4) 

The best fit of this to k{ijj) in Fig. [T] indeed requires e = 
1.8510“^^ erg/rad^ and is drawn as a line superimposed 
on the dots. The maximum angular frequency obtains for 
ka/2 — 7r/2 (i.e. when the phonon wavelength reaches its 
minimum, 2a) and is then 


cc;max = = 1.3110^^ rad/s. (5) 

ma^ 

The inflexion point ojc = cj/2 is at 0.66 10^^ rad/s, where g is 
minimum and V maximum, as confirmed by Fig. [T] The in¬ 
terpretation of the modeling results for transverse vibrations 
is therefore validated. The physical reason for the variation 
of V is that the degree of bending of the molecule depends 
on k and, therefore, on u. 

The same procedure can be applied to the longitudinal 
vibrations shown in Fig. [3] Here the stretching potential in¬ 
volved is quadratic in the variable a: AU — jSA . However, 
the fit of k{uj) now requires 


2 ^ . 2 / . 
muj = 4y5sm ( —), 


( 6 ) 


as in the classical treatment of solids, with /3 = 1.1210® 
erg/cm^ (617 kcal/mol/A ^) and Ct;max = 4.810^^ rad/s in 
the present instance. While the near constancy of V at low 
frequencies naturally derives from this equation, the near 
constancy of g is intriguing. The reason for it will become 
clear below. 

Note that, in both the transverse and longitudinal cases, 
the power index of k is even, indicating the existence of 
two waves traveling in opposite directions for each resonant 
frequency. 

In dealing with more complex molecular structures, it is 
no longer possible, by simple monitoring of each vibration on 
the computer screen, to distinguish phonons belonging to the 
same family, determine their wave vector, /c, and apply the 
same procedure. The sole quantity that can still be clearly 
defined and computed is the DOM. This difficulty can be 
circumvented as described next. 


3 GENERAL RELATIONS FOR 
N-DIMENSIONAL MOLECULES 

Consider first a periodic structure in n dimensions, with size 
L in every dimension. Let ki be the wave vector of vibra¬ 
tion i. According to the Born-von Karman (cyclic) boundary 
conditions, the structure can only sustain waves such that 
its length is an integral number of half-wavelengths. Then, 
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the number of possible vibrations, or states of the system, 
with k < /cat is 


N = 


2kN ^n 
27r/L^ ’ 


(7) 


where the factor 2 in the numerator accounts for the fact 
that k can be positive as well as negative. Then, the density 
of modes around kN (in /c-space) is 


. dN 2^kl- 

w[kN) = -j-j — — 


and the DOM is 


dkN (2 7r/L)^’ 

w{kN) 


= 


duj/dk' 


Hence, 


g{oj)(ko = Ck^~^dk ; C = n(-)". 


( 8 ) 


(9) 


( 10 ) 


Let G{cj) = f g{u})dcj, which is also the phonon rank when 
the phonons are sorted in order of increasing frequency. 
Then 


and 


O 


,,, , dio nLG^-^!'^ 
= 177: = 


dk 7T g 


(11) 


( 12 ) 


Thus, k and, more importantly, V, can be derived directly 
from the DOM, which, for large, complex molecules, is the 
only readily computable quantity. Note that g and G are pro¬ 
portional to the total number of modes N, which is known 
to be nearly equal to 3 times the number of atoms (for large 
molecules). This, in turn is proportional to the “volume” 
of the molecule (for molecules of the same structure), i.e. 
the product of all the dimensions,Li...Ln. Therefore, in the 
general case where all dimensions are not equal in size, L 
in eq. 12 should be equated with the geometric mean of the 
dimensions, (Li...Ln)^. In any case, the dependencies of g 
and G on L cancel out the factor L in Eq. 12, so the phonon 
velocity is in fact independent of size, as it should. On the 
other hand, the longest wavelength sustained by the struc¬ 
ture is not determined by L but by the longest dimension of 
the structure. 

This treatment is only valid in spectral ranges where no 
overlapping of vibrational modes of different types (trans¬ 
verse/longitudinal) occur, which is only the case for the 
lowest-frequency transverse phonons. In that case, the re¬ 
sults of the theoretical treatment for the speed of transverse 
waves agree with those of chemical modeling to better than 
a factor 2. This difference itself can be traced back to the ab¬ 
sence of the correcting factor / (Sec. 2.1) in the calculations, 
because of the assumed Born-Karman conditions. 

Now, since we are only interested in the lowest frequen¬ 
cies, we are entitled to set, as suggested by the log-log graph 
in Fig. 2, and as an approximation. 


g{(jj) (X V{(jj)(xuj^. (13) 

It follows that 

F(w) (14) 


^ 10 - 



E 



1 


10 


,(A°) 


50 


Figure 4. Dots: minimum wavenumber as a function of 
the largest molecular dimension. Line: power law fit, z/ = 
2U8/p-^t° '^^ cm-i = 6.4104/Li y°-25 GHz. The linear car- 
bon chains which served as a basis for the theoretical analysis 
(Sec. 3 and 4) were not included in this graph as they do not 
seem to be abundant in space. However, they too fall well within 
the upper and lower limits of the power-law index in Eq. 18. 


and, finally. 


m 1 

p -\ -1-= 1. 

n n 


(15) 


Thus, in 3D (n = 3) homogeneous media, the velocity 
is constant (as are the density and elastic constants), so 
p = 0 and, hence, m = 2, as first shown by Debye. On 
the other hand, for the ID chains considered above n = 1 
so m = —p, as found for both transverse and longitudinal 
vibrations, while both are 0 for the longitudinal ones. There 
is no apparent reason why n should be restricted to integers. 

In fact, both p and m are determined not only by the 
dimension n, but also by the dominant force governing the 
mechanics of the structure. The latter translates into the 
power index, s, in the relation k (x for low frequencies, as 
in Eq. 4 (bending force) and 6 (stretching force), which cover 
the more usual cases for common molecules. In these cases, 
s = 0.5 and 1, respectively. An immediate consequence is 


p = 1 — s, so m = ns — 1. 


(16) 


The fundamental, independent parameters of the problem 
are n and s. While the DOM depends on both, the veloc¬ 
ity power index does not depend on n, which may take any 
value between 1 and 3. As a consequence, while m can have a 
wide range of values, p is essentially constrained to two val¬ 
ues, 0.5 and 0, corresponding to transverse and longitudinal 
vibrations respectively. 

Non-integer values of n occur when, for instance, the 
structure consists in chains of catacondensed single carbon 
rings; n then lies between 1 and 2, because, strictly, the 
structure is neither ID nor 2D. In such cases, the constraint 
on p helps lifting the ambiguity in n. 

Some of these results are illustrated next. 
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4 MINIMUM AND MAXIMUM PHONON 
FREQUENCIES 

The minimum (or fundamental) phonon frequency, z/min = 
z^(l), is that for which k = Tvf/Lmoi, where / is the cor¬ 
rection factor for boundary effects (Sec. 2). Now, if k oc 
(Sec. 3), then it follows that — uj/sk{ijj). Therefore, 


Ki) 


sfV{\) 

2cZ/rxiol 


(17) 


The corresponding mode is always that of a transverse 
phonon, because the stress governing longitudinal waves is 
stronger than that which governs transverse waves, making 
their velocity higher. If the velocity were independent of 
frequency (s = 1), z/min would simply scale like the inverse 
of the molecular size, but this condition is far from common, 
except for 3D structures; as indicated above, for TA modes 
in ID and 2D, the common value of s is 0.5. 

In order to study relation 17 for more complex molecules 
than the linear ones in Sec. 2, 15 molecules were built, with 
backbones made up of 5- and 6-membered carbon rings, 16 
to 42 A in major dimension (see Appendix). Figure U col¬ 
lects the minimum wave numbers delivered by their normal 
mode analysis, as a function of molecular size. Although the 
structures are quite different and most of them are not plane, 
but bent and warped, a clear trend is apparent. The linear 
fit shown is 


= (18) 

in cm“^ and A , which is good to a factor 2. For v in GHz, 
the pre-factor would be 6.4 10^. 

On the other hand, inserting the same data (Lmoi, 
directly into Eq. 17, and plotting Umin against z/min, one 
obtains Fig.[5l Superimposed upon the dots is a linear fit to 
this plot: 

T 7- _ o CT 1 0.55it0.06 /-i pv\ 

Umin = 3.510 Z/^in 5 (19) 

whose exponent is close to the expected one, in agreement 
with the finding in Sec. 3 that p = 1 — s is constrained 
to 0.5 for transverse waves. The “average” velocity is close 
to the velocity found for linear molecules, indicating that 
the dominant force in both ID and 2D is that which resists 
bending, and that boundary effects are limited. Application 
of Eq. 12 to individual molecules confirms this result, as 
illustrated by examples in the following Section. 

If p = 0.5 is assumed to be the exact value, then Eq. 
17 leads to z/min cx: with which Eq. 18 agrees within 

error bars. A number of causes may contribute to the dis¬ 
persion of data points in Eig. [3]and[5l a) Lmo\ is not clearly 
defined a priori^ unless the molecule is definitely a straight 
line or a plane surface; b) z/min is uncertain because of larger 
code errors at low frequencies and of mode coupling with 
fingerprint or longitudinal vibrations; c) boundary effects 
and uncertainty on the correction factor /; d) the velocity 
at a given frequency depends on mass (through m), density 
(through a) and bonding force strength (through e or /3), all 
of which depend on the particular structure. 

The lack of computing power of presently available 
chemical modeling codes makes it difficult to explore larger 
structures. However, a conservative estimate may be ob¬ 
tained using Eq. 18 or Eig. [U Thus, in order to contribute 



Figure 5. Dots: Transverse acoustic sound velocity derived from 
eq. 17, for the lowest phonon frequencies of the 15 carbon-based 
molecules of Fig. 0] . Line: power law fit, V = 3.5 10"^ 
cm/s. 


absorption/emission around 100 GHz and beyond, a struc¬ 
ture has to be larger than about 40 A . 

For a periodic chain of N atoms, the maximum phonon 
frequency corresponds to a vibration wavelength about equal 
to the atomic periodicity, a. For a TA phonon on a carbon 
chain, this is of order 1000 cm“^, as in Fig. [U which may 
be extended to carbon structures in general. In between the 
minimum and maximum frequencies, there are about 3N 
modes. Any slight change in shape or composition of the 
structure entails a slight redistribution of the whole spec¬ 
trum. If the line of sight of a spectrometer includes enough 
molecules of the same structural type, then the spectrum 
will take the appearance of a continuum from about 1000 
cm“^ downwards. As suggested by Papoular (|2014ap . this 
might apply to space, where huge quantities of molecules 
form around stars, novae and supernovae in a haphazard 
way but in the same local environment. And, indeed, when 
observing galaxies or planetary nebulae one obtains contin¬ 
uous emission spectra starting from about 1000 cm“^ and 
extending down to 1 cm“^ or less: see Smith et al. ()2007D . 
Gompiegne et al. ()201ip . Fischer at al. (|2014p . Zhang and 
Kwok ()2011l) . 


5 SOME PARTICULAR ATOMIC 
STRUCTURES 

To illustrate finite 2D structures, three compact benzene 
clusters, 17x17, 13x17 and 7x34 A ^ in size, respectively, 
were built, optimized and analyzed for normal modes. The 
lower panel of Fig. [6] displays the DOM, g{i/) of the 17x17 
A ^ structure, computed from its spectrum delivered by 
the simulation code. Here, the spectral range of pure TA 
phonons is limited downwards by the structure finiteness 
and and upwards by boundary effects (edge hydrogens). 
Beyond that, the overlap of phonons of different types, as 
well as fingerprint modes, make interpretations hazardous. 
A straight line g = 0.22 cm is superimposed upon the DOM 
data, implying m = 0, in agreement with Eq. 16 with n = 2 
and s = 0.5. 
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Figure 6. Benzene cluster 17x 17 A 2 in dimensions. Bottom: 
DOM, delivered by the chemical model; line: power law fit 

to the phononic tail, g = 0.22 cm. Top: same conventions for 
V{uj)] line: power law fit to the phononic tail, V = 3.310^ 
cm/s. 



Figure 7. Chain of 11 anthracenes and phenanthrenes (126 C, 66 
H). Bottom: DOM, g(u)^ delivered by the chemical model; line: 
power law fit to the phononic tail, g = cm. Top: same 

conventions for line: power law fit to the phononic tail, 

V = S.llO^z/O-^ cm/s. 


The upper panel displays the corresponding velocities, 
as deduced from g(iy) and eq. 12. The straight line represents 
a power law: V = 3.3 10^ cm/s, so p = 0.5. Note that m 
and p exactly satisfy eq. 15 for n = 2. This is not the case of 
the DOM and velocities of the other two benzene clusters, 
for their 2 dimensions are not exactly equal. If the average 
DOM of the 3 structures is considered, its low-frequency tail 
indicates m = —0.1. Inserting this value in Eq. 16, together 
with s = 0.5, one obtains n = 1.8, a measure of the excur¬ 
sion from the proper 2D case. This excursion increases, as 
expected, with the aspect ratio of the structure. Thus, for a 
chain of 11 anthracenes and phenanthrenes (126 C, 66 H), 
one finds m = —0.3, so, with s = 0.5, we get n = 1.4 ( Fig. 
0 . 

The transverse vibrations of 2D molecules is reminis¬ 
cent of plate and drum vibrations. However, it obeys, in 
fact, different mechanics because, whatever the thinness of 
a solid state sheet, the shear and compression will remain 
the dominant stresses. 

As an example of unexpected findings, consider a chain 
of six “trios” connected in a row, one to another, by a sin¬ 
gle C-C bond. Each trio is made of one 5-membered carbon 
ring squeezed between 2 6-membered rings, the free vertex 
of the pentagon being occupied by a Sulfur atom. Each C-C 
connection is subject to bending and torsion. Despite the 
planar character of each trio, relation 15 is only satisfied if 
the system is treated as one-dimensional, using the proce¬ 
dure of Sec. 3. Figure [8] displays the results. In each graph, 
an average power-law line is superimposed upon the dots. 
The power-law indexes are -0.5 and 0.5, respectively, as for 
transverse phonons in linear chains, but the velocities are 
about a factor 0.7 lower than those of phonons in linear 
chains, at the same frequencies. 



V (cm’^) 

Figure 8. Chain of six “trios” connected in a row, one to another, 
by a single C-C bond; each trio is made of one 5-membered carbon 
ring squeezed between 2 6-membered rings. Bottom: DOM, g(n)^ 
delivered by the chemical model (dots); line: power law fit to the 
phononic tail, g = 0.2 cm. Top: same conventions for V{u)] 

line: power law fit to the phononic tail, V = 2.2 10^ cm/s. 


6 EXTINCTION CROSS-SECTION OF 
PHONONS 

In order to assess the spectroscopic detectability of these 
phononic vibrations, it is necessary to evaluate the inten¬ 
sity of the attending electromagnetic signal, if any. This de¬ 
pends, of course on the particulars of the molecules at hand. 
However, some general properties can be identified. The fol¬ 
lowing study is restricted to carbon-based molecules which 
are 1- or 2-dimensional in the broad sense (i.e. not necessar¬ 
ily straight or plane, which is the case of most molecules in 
space). 

After optimization of the structure, the modeling soft¬ 
ware that is used here performs a Normal Mode analysis 
which delivers all mode frequencies with the correspond¬ 
ing integrated absorption intensities, I(iy) in km/mole (1 
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Figure 9. The integrated absorption intensity for the transverse 
phonons in linear carbon chains of 20 (red dots) and 160 atoms 
(black +). 



Figure 10. The integrated absorption intensities for a benzene 
cluster 17x17 A ^ in dimensions. 


mole=6.024 10^^ molecules) . This is related to the absorp¬ 
tivity, a, of a gas of identical molecules by 

a{iy) = 10^ I ( 20 ) 

where a is in cm“^, C is the molecular density (mole/1) and 
Az/, the band width of the band at u, in cm“^. For a large 
assembly of molecules of the same family but slightly differ¬ 
ent from one another, so as to produce a quasi-continuous 
vibrational spectrum, it can be shown (Papoular 2014b) that 

a{iy) = 10 ^ < I{iy) > g{iy)C. ( 21 ) 


For astrophysical applications, it is often more convenient to 
consider the extinction cross-section per H atom, cr, which 
is related to a by 


aR) = /A 


«(0 

Vc ^T-dust 


/A < IR) > gR) 
6 IQi® Nc 


( 22 ) 


where a is in cm^, / is the fraction of carbon atoms locked in 
dust, A is the cosmic abundance of carbon relative to hydro¬ 
gen, Nc is the average number of C atoms in one dust par¬ 
ticle and Ud is the density of dust particle in space (cm“^). 

We now set out to estimate a for wavelengths beyond 
about 1000 /xm. Consider, hrst, the intensities, 1. Figure [9] 
displays I{iy) for the shortest carbon chain studied in Fig(T] 
( with CH terminations), and for the longest chain in the 
present study: 160 C atoms (terminated by a couple of H 
atoms at each end). 

Two categories of transverse phonons can be distin¬ 
guished: a) those with / between about 0.01 and 0.1 km/mol; 
b) the “suppressed” ones, with I one or several order of mag¬ 
nitude weaker. The distinction and the gap in / are particu¬ 
larly clear for the shortest chain. A careful study of this case 
shows that the strongest (respectively the weakest) intensi¬ 
ties are associated with odd (respectively, even) numbers 
of half-wavelengths along the chain, corresponding to con¬ 
structive (respectively, destructive) combinations of atomic 
excursions. This effect weakens as the chain length increases, 
progressively reducing the range of intensities from below. 

A third phonon category, the “enhanced” phonons, is 
discernible in the hgure, as the wave number increases, but is 
more obvious when the molecular geometry favors the edges 


relative to the backbone or when the edges involve more 
infrared-active functional groups like OH for instance. In 
ah these instances, the edge atoms dominate the molecular 
polarizability. 

Similar trends are observed in FiglTOl obtained from 
the modeling of the 17x17 benzene cluster of Fig[6l In 
this case, however, the absorption intensities of the low- 
frequency phonons are an order of magnitude higher and 
more and more high intensities occur as the frequency in¬ 
creases. Both effects are due to the increased number of 
atoms and to the presence of a larger relative number of 
edge hydrogen atoms. 

From the consideration of these, and the other struc¬ 
tures modeled in this work, and in view of the computa¬ 
tional difficulty of increasing sizes, it may tentatively be 
concluded that the intensities of the lowest-frequency (trans¬ 
verse) phonons, for carbon based molecules in general, tend 
to hock around a constant value as the frequency decreases. 
This value lies between 0.01 and 0.1 km/mol for molecules 
with n = 1, and 0.1 to 1 for structures with n — 2. 

As for the density of modes, g(y)^ which also appears 
in Eq. 22, it can be derived using the results of Sec. 3 and 
4. From Eq. 12 


and, from Eq. 13, 


(23) 


g{uy)=g{co,)Nr 


flL / ^ \m 
'kV{uJi) UOl 


(24) 


where cji is the lowest phonon angular frequency. Now, in 
the most common case of n = 2, then m = 0 and g does not 
depend on frequency. Erom Sec. 4, we have, approximately. 


1^1 


2148 

~L^ 


and V{iyi) 




1.6210® 

Z 


(25) 


so whith g in cm units and L in A , 


g{iy)=g{iyi) = 7A10*L\ (26) 


in which, we note in passing, the particular case of Eig. 6 
hts fairly well, as it should. Equation 22 also features Nc, 
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which is found empirically to be (Ll2)‘^, so 

^ = 310“’^ cm. (27) 

Nc 

Inserting this result into Eq. 22, together with A = 1/3000, 
and < I >= 0.1 (see comments above in this section), one 
finally gets 


(jc —1.710 / cm^/Hatom. (28) 

If one takes / = 0.15, as did Papoular (|2014bD to model the 
high Galactic latitude diffuse emission, then a ^ 2.610“^^ 
cm^/H atom, independent of frequency and slightly weaker 
than the value found in that work for the longest wave¬ 
length, A = 1000 fim This confirms the preliminary conclu¬ 
sion drawn in Papoular (|2014ap . to the effect that the ex¬ 
tinction cross-section of carbon-rich dust particles levels off 
beyond a few hundred micrometers. It must be recalled that 
this extinction estimate is only valid to as long a wavelength 
as allowed by the maximum size of the particles present in 
the dust (Eq. 18). Towards shorter wavelengths, extinction 
progressively increases as functional group modes set in with 
increasing density and intensity, and the extinction settles 
into a logarithmic power law with index [3 between 1 and 2, 
depending on environment. 


7 EMISSION SPECTRA 

Knowing the integrated intensity spectrum of a molecule, it 
becomes possible to compute its emission spectrum if the 
excitation process is also known. The simplest case is that 
of thermal equilibrium with the environment at a given tem¬ 
perature. The emission spectrum is then obtained by multi¬ 
plying the Black Body spectrum with the absorption cross- 
section (Eq. 22). An example of this procedure is given be¬ 
low. However, for this to be valid, the excitation rate must be 
strong enough to reach statistical equilibrium between emis¬ 
sion and excitation processes. Eor low intensities of ambient 
radiation (G=Go), this implies that the absorption cross- 
section of the molecule is large enough (size larger than 0.1- 
1 /im). Otherwise, before the molecule radiatively relaxes, 
it has time to settle in a state of internal statistical equi¬ 
librium, wherein the deposited excitation energy is continu¬ 
ously flowing between vibrational modes, but the energy in 
each mode remains constant on average. Roughly speaking, 
its emission spectrum will cover the range between 20 /xm 
and Amacc, the maximum wavelength in the phonon spec¬ 
trum: this range will often turn out to be much wider than 
that of a black or gray body. 

To illustrate this case, we select 3 organic molecules 
built up of 1, 3 and 6 “trios” (see the Appendix) arranged 
in a row so that their lengths turn out to be 9, 25 and 50 A 
, and their Xmax to be 369, 1275 and 5097, respectively, as 
shown in Eig. [TT] which displays their integrated absorption 
spectra. They are meant to typify groups of dust particles of 
similar structure and sizes. The peak of integrated absorp¬ 
tion which protrudes between 20 and 60 /xm is a fingerprint 
of the OH groups attached to each trio, which are expected 
to be common in the ISM . 



Figure 11. The integrated band intensities of the 3 selected 
molecules: 1 (red dots), 3 (olive triangles) and 6 “trios” (blue 
stars), respectively 9, 25,51 A in length. The abrupt end of each 
computed spectrum is indicated by vertical dashes. The maxi¬ 
mum mode wavelength is 5097 /xm. 


7.1 The computational procedure 

Now, it was shown (Papoular 2012) that the energy density 
ultimately radiated by an isolated molecule after it has been 
excited by radiation or collision that deposited an energy 
Eexc in it, and after it has promptly relaxed into internal 
statistical equilibrium, is 


dE*{iy) = Eexc 


e{iy)iy^dm{iy)^du 
T>e(iy)iy^dm(iyydiy ’ 


(29) 


where e{ip) is the fraction of energy staying on average in 
mode IP and m is the total electric dipole moment of the 
molecule. It should be remembered that the latter (in the 
ground state of the molecule) is linked to the integrated 
band intensity, I (km/mol), by 

I{ip) = 2.6ipm{ip)^ , (30) 


stressing the intimate connection between emission and ab¬ 
sorption. The modeling package used here provides for com¬ 
puting both distributions e{i') and These are obtained 

in 3 successive runs of molecular dynamics: 

a) After optimizing the ground state of the selected 
molecule, the latter is given an energy Eexc and left in vac¬ 
uum to relax into its internal statistical equilibrium, which 
is reached within a few tens of picoseconds. 

b) A second run is then launched, during which the fluc¬ 
tuations of the total kinetic energy in time are monitored. 
The Eourier transform of this delivers half the steady state 
kinetic energy distribution (e(z/)) among the modes, the 
other half being in the form of potential energy. This run 
must be at least twice as long as the period of the lowest 
frequency phonon. 

c) A third run is launched to monitor the temporal fluc¬ 
tuations of the total dipole moment. The Eourier transform 
of these delivers the m(iy) spectral distribution. 
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Figure 12. Square dots: Diffuse High Galactic Latitude SED, 
adapted from Compiegne et al. (I2011D . Circular dots: example 
of measured AME, adapted from Abergel et al. il2nI3t . Spectra 
computed here for the three molecules made up of: 1 (red line), 
3 (olive dotted line) and 6 (blue dash-dotted line) “trios” in a 
row, assuming G=Go and w = 1 (see text) in each case. Note the 
horizontal trend beyond the fingerprint range and the similarity, 
in relative values, of these emission spectra with the correspond¬ 
ing absorption spectra of FigElThe abrupt end of each computed 
spectrum is indicated by vertical dashes. It would take a molecule 
of 15 “trios” in a row (390 atoms and 125 A in length) to cover 
the whole range of observed anomalous emission.. 

7.2 Results 

Assume the rate of excitation events corresponds to Go, or 
10®cm“^s“^, and the total absorption c-s (cross-section) of 
all interstellar dust is 10“^^ cm^ per H atom (Spitzer 1978). 
The total rate of excitations is therefore 10~^^w s“^, where 
w is the fraction of dust cross-section contributed by a se¬ 
lected particle type. The emission power is obtained by mul¬ 
tiplying this with dE^*(z^] (Eq. 29); it is shown in Fig. [12] for 
the 3 selected structures, assuming w — 1 and 10^° H atoms 
per cm^. The close similarity, in relative values, of these 
emission spectra with the absorption spectra (Fig.fTTl) is re¬ 
markable and suggests that the latter are valuable guides 
in the study of emission. Another remarkable feature of the 
spectra is their flatness beyond the fingerprint region (>100 
tiro). 

It has not proved possible to study longer molecules, 
as presently available modeling codes cannot handle parti¬ 
cles larger than a few hundred atoms. However, the spectral 
extension to longer wavelengths clearly increases with the 
length of the molecule as suggested by Sec. 4 and illustrated 
in Fig. [TTl Extrapolation of the present trend of Xmax with 
molecular length indicates that merely increasing the latter 
to 130 A would suffice to cover the whole AME (Anomalous 
Microwave Emission; see Abergel et al. 2014, Genova-Santos 
et al. 2015) spectral range, to 10 GHz. 

Figure [12] also displays one version of the observed 
emission spectrum from the diffuse high latitude Galactic 
medium (adapted from Gompiegne et al. 2011 and Abergel 
et al. 2014). Gomparison with the spectra computed here 
suggests that, with similar carbon-rich molecules, assuming 
G=Go, only a small fraction, re, of the total dust is needed 
for their emission to reach the levels observed in the high lat¬ 
itude Galactic medium, between 20 and 60/xm, and beyond 



Figure 13. Emission of a 6-trio molecule under thermal equilib¬ 
rium at 20 K (red dashes), or isolated in vacuum after thermal- 
ization at 20 K (black line). 


the 18-K gray-body spectrum (Gompiegne 2012). Similar 
molecules include, of course, graphitic particles like PAHs 
and graphite. 

The evolution of the emission spectrum as a function of 
the molecular length illustrated in Fig.[T2]can be understood 
in the light of the phonon physics developed in Sec. 3 and 
4, as shown in Appendix B. This also provides the grounds 
for extrapolation towards larger molecules and longer wave¬ 
lengths. 

Figure [13] compares, for the longest structure (six trios 
in a row), emissions at low and high (thermal equilibrium) 
rates of excitation, or, correspondingly, small and large par¬ 
ticles. This highlights the difference in spectral extensions 
of emission for the same molecule. 

The present results are potentially of use in the on¬ 
going discussion of the carriers of AME and the degree of 
polarization of their emissions as compared with the GMB. 
This problem recently gained traction in the wake of the 
measurement, by the BIGEP2 telescope, of B-mode emission 
polarization in the spectral range of the cosmic microwave 
background (GMB), which could be a signature of early cos¬ 
mic inflation (see Ade et al. 2014). The relevance, to this 
issue, of the molecules studied above depends first on their 
anisotropy. This may be illustrated by an example, viz. the 
6-trios molecule of Fig. [12] Our chemical modeling package 
finds the 3 principal axes of inertia to be, respectively, 23036, 
180748 and 192842 a.m.u. A which also define the coordi¬ 
nates x,y and z. The mechanical anisotropy is obvious. The 
permanent dipole moment is 4.34 Debye, with components 
1.98, -3.03 and 2.4. 

In order to assess the anisotropy of the vibrational emis¬ 
sion, one follows in time each of the 3 components of the 
molecular dipole moment, then takes the Fourier transform, 
Ii(y) (i:x,y,z) of its auto-correlation function. This delivers 
3 curves similar and parallel to the blue curve in Fig. 1141 of 
Appendix B. Now, after the molecule has been excited in a 
bath at 20 K (equivalent to a thermal spike), and left free 
in vacuum, it is found to rotate with an angular velocity of 
2.3 10“^ rad/ps (~0.37 GHz). This is well within the range 
considered by Lazarian and Draine (|2QQQ|) both for their 
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resonant relaxation and for Davis-Greenstein paramagnetic 
relaxation mechanisms. 

In the former mechanism, alignment in a magnetic field 
occurs through rotation relaxation because of spin-lattice in¬ 
teraction. This requires that a quantum of rotation energy 
may be absorbed by the lattice, i.e. that resonance may oc¬ 
cur between rotation and one of the lattice vibrations, for 
instance. For the 6-trios molecule considered here, this is 
not possible because the lowest vibrational frequency in the 
present example, is 60 GHz, much higher than the rota¬ 
tion frequency, 0.37 GHz. If a bigger molecule is considered 
or some mechanism operates to give this molecule a higher 
rotation energy, then spin-lattice relaxation might become 
effective. 

Anyhow, assume some alignment mechanism is at work 
and forces the molecular rotation vector to lie parallel with 
the Galactic magnetic field, which, in turn is supposed to 
be perpendicular to the sight line. In this picture, there¬ 
fore,one component of the dipole moment lies along the 
magnetic field while the other two are randomly oriented in 
the plane perpendicular to this direction and containing the 
sight line.Three polarization fractions can then be deduced, 
according to which dipole moment component is parallel to 
the magnetic field: 


Pi{v) =1 


Ijiy) + Ikiy) - h{v) 

Ij(v) + Ik(v) + liiv) 




(31) 


Plotting the 3 quantities as a function of frequency, it is 
found that they oscillate between 0.4 to 0.95 at the phonon 
frequencies and nearly zero in between. Gonsidering a fam¬ 
ily of similar molecules, differing only in length, and thus 
increasing the density of phononic states, one expects each 
polarization to retain high values throughout the spectrum. 
It must be stressed, however, that this is assuming perfect 
alignment and rigid structures. This whole exercise is only 
meant to show the potential of this type of molecules as 
polarized emitters in the AME spectral range. 


8 CONCLUSION 

Several neutral, carbon-based, hydrogen-terminated, 
energy-optimized molecules, made of chains and rings, 
were modeled by means of a semi-empirical algorithm. 
Their building blocks are simple enough that they may be 
synthesized in the laboratory and in space. 

Their vibrational spectra were obtained by Normal 
Mode analysis. The focus was set upon those vibrations that 
extend over the whole structure (phonons), setting aside fin¬ 
gerprints due to small, inevitable, functional groups (con¬ 
tributing to the high frequency end of the spectrum). Trans¬ 
verse and longitudinal acoustic vibrations of the simplest 
structures could be separated by monitoring the atomic mo¬ 
tions of each mode. The asymptotic behavior of their DOM 
and sound velocity toward low frequencies were found to 
obey power-laws in the frequency, with indexes m and p, re¬ 
spectively. For the transverse vibrations, these indexes differ 
from the solid state values, 2 and 0. It was shown that this 
is so because the force governing the motion is that which 
resists bending of G-G bonds. The intensity of this force was 
confirmed independently. Because this force is weaker than 


the stretching force of G-G bonds, the lowest frequencies in a 
phonon spectrum are always those of TA vibrations. Knowl¬ 
edge of the asymptotic TA phonon velocity, in particular, is 
essential in determining the lowest frequency to be expected 
in the spectrum of a given molecule. 

In the case of more complex structures, it is hardly pos¬ 
sible to separate phonons of different types. The only quan¬ 
tity readily available from modeling is the DOM. Section 3 
gives a general formula from which the sound velocity can 
be derived, given the DOM. 

Application of this formula requires the knowledge of 
the dimensionality, n, of the structure. For common struc¬ 
tures, warped and finite, this notion is not as straightforward 
as it appears to be: n is not always an integer, 1 or 2, but 
rather lies in between. It was shown that the ambiguity can 
be lifted by using a general relation between m, n and p. 

These results were illustrated by the study of several 
structures of different types (but not crystalline). It was 
found that the TA sound velocity always scales like 
so the lowest frequency in a spectrum scales like where 
Lmoi is the longest straight dimension of the structure. The 
extinction cross-section of carbon-rich particles was shown 
to level off progressively beyond about 100 pm (<3000 GHz), 
and was estimated at 3 10“^^ cm^/H atom, in order of mag¬ 
nitude, in this spectral range. 

The results obtained above can readily be used to de¬ 
duce emission spectra, assuming the ambient radiation den¬ 
sity is high enough (or the dust particle large enough) that 
it reaches thermal equilibrium before re-radiating the en¬ 
ergy gained. In the opposite case, the emission spectrum 
can be deduced from the monitoring of the temporal fluctu¬ 
ations of the total electric dipole moment after the molecule 
has been excited . From spectra computed in this way for 
molecules up to 50 A in length, it is concluded that rea¬ 
sonable amounts of G atoms locked in similar structures are 
sufficient to emit detectable radiation as far as 4000 pm at 
least. Presently available modeling codes can hardly handle 
larger particles, but the theoretical analysis indicates that 
molecules 120 A in length are needed to cover the whole 
spectral range thought to be dominated by dust emission, 
before free-free bremsstrahlung and synchrotron emission 
prevail. 

The primary application of this work is to the exten¬ 
sion of the interstellar extinction curve and the variation 
of its slope. The relation between molecular size and low¬ 
est phonon frequency should help define the upper size limit 
of dust particles and perhaps their size distribution as well. 
The present results also shed new light on the vibrations of 
molecular structures that are not commonly studied. It was 
also shown that hydroxyl-bearing G-rich molecules can pro¬ 
vide the emission observed between 20 and 60 pm, which is 
usually assigned to previously undefined very small grains 
(VSG; see Desert et al. 1990). 

One could also envision, in the future, an application 
to the elucidation of the sources and properties of AME. 
The challenge is to separate the contributions of GMB and 
foreground dust to the measured B-mode polarization (see 
Adam et al. 2014). Another challenge is to determine the 
cause(s) of emission polarization in the spectral range of 
interest. The present work indicates that large carbon-rich 
molecules may contribute significant polarized phonon emis¬ 
sion in that range. Because of the intricacy and intertwin- 
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ing of CMB and dust emissions, an accurate quantitative 
knowledge of dust polarization is required. This is certainly 
an interesting subject of further investigation. 
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9 APPENDIX A 

Table 1 provides the atomic composition, and the values of 
Lmoi and z^min of thc molecules used in building Figjd] and 
[5l in order of increasing ^min • 

Here is a brief description of their structures: 

1) 6 “trios” connected by single C-C bonds; each trio is 
made of a 5-membered ring squeezed between 2 6-membered 
ring; each has 1 edge hydroxyl and one -CH20H; 

2) simulation of a macromolecule of kerogene; 

3) same as (1) but O atoms excluded; 

4) same as (1), except: no hydroxyls, only two -CH20H; 

5) 5 “trios” connected by single C-C bonds; each trio is 
made of a 5-membered ring squeezed between 2 6-membered 
ring; each has 1 edge hydroxyl and one -CH20H; 

6) 4 “trios” connected by single C-C bonds; each trio is 
made of a 5-membered ring squeezed between 2 6-membered 
ring; each has 1 edge hydroxyl and one -CH20H; 

7) linear chain of 15 benzenic rings; 


[ht] 


Table 1. The structures represented in Fig. [4| and [5] 


No 

Nat 

G 

H 

O 

S 

Tmol 

^min 

1 

146 

78 

50 

12 

6 

42 

1.96 

2 

493 

254 

185 

40 

9 

35 

2.34 

3 

119 

73 

40 

0 

6 

42 

2.65 

4 

130 

76 

46 

2 

6 

42 

2.67 

5 

122 

65 

42 

10 

5 

32 

3.76 

6 

98 

52 

34 

8 

4 

48 

1.27 

7 

96 

62 

34 

0 

0 

38 

4.5 

8 

142 

103 

38 

1 

0 

37 

4.79 

9 

101 

51 

39 

8 

2 

25 

6.15 

10 

62 

37 

22 

0 

3 

20 

7.4 

11 

59 

10 

22 

0 

0 

23 

8.4 

12 

292 

110 

182 

0 

0 

16 

11.4 

13 

152 

120 

32 

0 

0 

17 

13 

14 

74 

49 

24 

1 

0 

20 

15.2 

15 

114 

90 

24 

0 

0 

17 

19.8 

The units of Lmol cind z/min 

L are 

A 

and cm“ 

respectively. 


8) linear chain of 4 coronenes; 

9) 1 pyrene, 1 trio, 1 benzenic ring; 1 pentagon; 

10) 3 trios, 1 edge CHS; 

11) 3 trios 

12) hydrogenated amorphous carbon (HAC); 

13) square cluster of benzenic rings (17x17 A ^); 

14) 2 coronenes, 1 edge -CH20H; 

15) cluster of benzenic rings (13x17 A ^). 


10 APPENDIX B 

We seek here to understand the physical origins of the slope 
and ordinate at the longer wavelengths in the spectra in Fig. 
M We also seek to predict the emission power of molecules 
that would be large enough to contribute to the observed 
signals at the longest wavelengths of the AME spectrum. 
Sections 3 and 4 are helpful for this purpose, and the same 
symbol conventions will be used. 

First, consider Fig. M which displays some of the inter¬ 
mediate results that were used in drawing Fig ll3l It suggests 
that, in the range 1-100 /xm, where the spectra are free of 
fingerprints features, the energy and electric dipole moment 
approximately follow power-law trends. Let the correspond¬ 
ing indexes be E and M, which are found to be -1 and -3. E, 
which characterizes the energy distribution among modes in 
statistical equilibrium in vacuum, after excitation, is a re¬ 
sult of the coupling between modes, and must be considered, 
here, as empirically given for it depends on the particular 
couplings and non-linearities of the molecule under consid- 
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Figure 14. Spectra of energy (red line) and dipole moment 
squared (blue dashes) for the six-trios molecule, thermalized in 
a bath at 20 K, then isolated in vacuum and monitored; arbitrary 
units. The corresponding power indexes are shown. 

eration. M, however, is linked to E by the relations obtained 
in Sec. 3. 

In transverse mode No z of a linear chain, the transverse 
atomic displacement is of the form yi = Yism{kira), where 
r is the rank of an atom along the chain. Then, from Eq. 1, 
the energy carried by the mode is 

Ei oc {kia/2)ErSiiP(kira) oc~ LY^ (kiaj'E)^ , (32) 

where L is the chain length and the sinus was approximated 
by its argument, as we are dealing with long wavelengths. 
Based on Sec. 3, ki oc (z/L); from Sec. 4, z oc . Now, 

the element of dipole moment squared, dmf, scales like Y^. 
Figure [m therefore suggests that we can put Y^ — Ah>P. 
Hence, 



(33) 


From Fig[T4l this scales like so F = 2 + M, which is 
satisfied by the indexes obtained empirically in that figure. 
Finally, Eq. 29 gives 


Z/dF*(z/) oc Eexcl^h'^l^^LAlJ^ oc Eexc , (34) 


independent of the frequency. That the emission spectrum is 
flat beyond the fingerprint region is very nearly confirmed, 
on average, by the modeling curves of Fig. 1121 These various 
instances of relatively good agreement between analytical 
and modeling predictions provide grounds for extrapolating 
the present results to larger structures that cannot be han¬ 
dled by present modeling codes. 




